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Abstract 

We establish symmetrization results for the solutions of the linear fractional dif- 
fusion equation dtu + (— A)'^/^u = / and its elliptic counterpart hv + (— A)°'/^z; — /, 
h > Q, using the concept of comparison of concentrations. The results extend to 
the nonlinear version, dtu + (— A)°'/^yl(u) = /, but only when A : M+ — >■ K_|_ is a 
concave function. In the elliptic case, complete symmetrization results are proved for 
B{v) + {-Ay/^v = / when B{v) is a convex nonneeative function for v > with 
-6(0) — 0, and partial results when B is concave. Remarkable counterexamples are 
constructed for the parabolic equation when A is convex, resp. for the elliptic equa- 
tion when B is concave. Such counterexamples do not exist in the standard diffusion 
case a = 2. 

1 Introduction 

The techniques of symmetrization are a very popular tool of obtaining a priori estimates 
for the solutions of different partial differential equations, notably those of elliptic and 
parabolic type. Symmetrization techniques appear in classical works like [29^ I36j . The 
application of Schwarz symmetrization to obtaining a priori estimates for elliptic problems 
is described by Weinberger in [56j, see also ^33j. The standard elliptic result refers to the 
solutions of an equation of the form 

Lu = f, Ln = - ^ di{aijdju) , 
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posed in a bounded domain U C M^; the coefficients {uij} are assumed to be bounded, 
measurable and satisfy the usual ellipticity condition; finally, we take zero Dirichlet bound- 
ary conditions on the boundary d^l. The now classical analysis, introduced by Talenti [40|, 
leads to pointwise comparison between the symmetrized version of the actual solution of 
the problem u{x) and the radially symmetric solution f of some radially symmetric 
model problem which is posed in a ball with the same volume as il. Sharp a priori esti- 
mates for the solutions are then derived. We refer to the papers |43t Hj in the framework 
of linear operators, and |42^ I27j where comparison results are obtained when dealing with 
nonlinear elliptic operators of divergence type. There is a large literature in this topic 
with many interesting developments. 

When this technique is applied to parabolic equations, the general program of comparison 
with a model problem of radial type still works, but the result of pointwise comparison need 
not hold and has to be replaced by comparison of norms at every time t > 0. Actually, 
a more basic result, called comparison of concentrations is true, cf. Bandle [71 [H] where 
linear parabolic equations with smooth coefficients are discussed. Such results have been 
extended in works like [35] , [5] , |55] for weak solutions of linear parabolic problems with 
discontinuous coefficients. Relevant definitions about symmetrization, rearrangements 
and concentration are recalled in Section [2l 

Elliptic approach to nonlinear parabolic problems. An extension of the sym- 
metrization results to nonlinear parabolic equations of possibly degenerate type, more 
precisely of the porous medium type, was done by the first author in ^45j. The paper 
considers the evolution problem 

(1.1) dtu = AA{u), n(0) = no, 

where A a monotone increasing real functioii[^ and uq is a suitably given initial datum 
which is assumed to be integrable. For simplicity the problem was posed for x G M^, but 
bounded open sets can be used as spatial domains. The novel idea of the paper was to use 
the famous Crandall-Liggett Implicit Discretization theorem |19j to reduce the evolution 
problem to a sequence of nonlinear elliptic problems of the iterative form 

(1.2) -hAA{u{tk))+u{tk) = u{tk-i), k = l,2,---, 

where tk = kh, and /i > is the time step per iteration. Writing A{u) = v, the resulting 
chain of elliptic problems can be written in the common form 

(1.3) hLv + B{v) = f, B = A-\ 

General theory of these equations, cf. [12], ensures that the solution map: 

T : f^u = B{v) 

^More generally, A can be a maximal monotone graph, but that generality is of no concern for us here. 
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is a contraction in some Banach space, which happens to be Note that the constant 

/i > is not essential, it can be put to 1 by scahng. In that context, the symmetrization 
result can be split into two results: 

(i) the first one applies to rearranged right-hand sides and solutions. It says that if two 
r.h.s. functions /i, /2, are rearranged and satisfy a concentration comparison of the form 
fi '< f2, then the same applies to the solutions, in the form B{vi) ~< B{v2)^ 

(ii) The second result aims at comparing the solution v of equation ( |1.3[ ) with a non- 
rearranged function / with the solution v corresponding to f^, the radially decreasing 
rearrangement of /. We obtain that t; is a rearranged function and B{v'^) -< B{v), i.e., 
B(v) is less concentrated than B[v). 

This precise pair of comparison results can be combined to obtain similar results along the 
whole chain of iterations u{tk) of the evolution process, if discretized as indicated above. 
This allows in turn to conclude the symmetrization theorems (concentration comparison 



and comparison of U' norms) for the evolution problem (1.1). This approach has had 
a large expansion in the past decades, cf. |17] and references. There is no difficulty in 
considering equations with a right-hand side, like Ut = A^(n) + /, as long as f £ L^{Qt), 
Qrp = X (0, T). We also mention how time discretization and symmetrization tools can 
be combined together to get interesting comparison results for some types of parabolic 
equations with double nonlinearity, such as b{u)t = /S.A{u) + /, with special assumptions 
on 6, see L20j and [6j, and to equations with weights [37j. The technique also applies for 
p-Laplacian operators, cf. [i6] . 



Equations with fractional operators. The study of ehiptic and parabolic equations 
involving nonlocal operators, usually of fractional type, is currently the subject of great 
attention. In that sense, it is quite natural to investigate how to apply symmetrization 
techniques to the elliptic equations like 

(1.4) {-^r/\ = /, 

where the standard Laplacian operator A is replaced by one of the fractional Laplacian 
operators (— A)"^/^, < o" < 2, as defined in [32l|39]. The study of this question has been 
successfully implemented by the second author and Di Blasio in a recent paper |21j by 
an interesting technique of analysis of Steiner symmetrization of an equivalent extended 
problem, based on the extension technique used by Caffarelli and Silvestre for the defini- 
tion of cr-Laplacian operator, [TCJ. For previous uses of Steiner symmetrization in standard 
elliptic problems see [3]. The results of [21] include a comparison of concentrations, in 
the form -< v, that parallels the result that holds in the standard Laplacian case; 
note however that no pointwise comparison is obtained, so the result looks a bit like the 
parabolic results of the standard theory as mentioned above. 

In the present paper we are interested in considering the application of such symmetriza- 
tion techniques to linear or nonlinear parabolic equations with similar fractional Laplacian 

^See the definitions of the order relation -< in Section [2] 
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operators. To be specific, we will focus on the equations of the form 



(1.5) dtu + {-Ay/^A{u) = f, 0<a<2. 

Following the theory for the standard Laplacian just sketched, we want to consider as 
nonlinearity A an increasing real function such that ^(0) = 0, and we may accept some 
other regularity conditions as needed, like A smooth with A'{u) > for all u > 0. The 
problem is posed in the whole space M^. We want to pay special attention to the form 
A{u) = u"^ with m > 0; the equation is then called the Fractional Heat Equation (FHE) 
when m = 1, the Fractional Porous Medium Equation (FPME) if m > 1, and the Frac- 
tional Fast Diffusion Equation (FFDE) if m < 1. Let us recall that the linear equation 
dtu + {-Ay/^u = is a model of so-called anomalous diffusion, a much studied topic in 
physics, see for instance [T| IHOl [Ml I54^ [57] and the references therein. The interest in these 
operators has a long history in Probability since the fractional Laplacian operators of the 
form (— A)*^/^ are infinitesimal generators of stable Levy processes, see [21 [111 HI]. 

For A(u) = n"* and general m > we obtain a nonlinear diffusion model; the theory of 
existence of weak solutions for the initial value problem has been addressed by the first 
author and collaborators in |22l 1231 l5T]. and the main properties have been obtained. In 
particular, if we take initial data in L^, then an L^-contraction semigroup is generated, and 
the Crandall-Liggett discretization theorem applies. Extension of this results to general 
smooth A is done in [25j. 

The application of the method of implicit time discretization leads to the nonlinear 
equation of elliptic type 

(1.6) h{-Ay/^v + B{v) = f 

posed again in the whole space or in an open subdomain C with zero Dirichlet 
boundary conditions; /i > is a non- important constant, and the nonlinearity B is the 
inverse function to the monotone function A that appears in the parabolic equation. 

Organization of the paper and main results. Section [2] contains the preliminaries 
about symmetrization and mass concentration that we will need. 

• As a first step of our analysis, we address in Section [3] the issue of comparison of 
concentrations for rearranged functions, more precisely how to compare the rearrangement 
of the solution of an elliptic problem with data / with the solution of a radial problem 
with data rearrangement of /. The technique used in does not work for the 



modified equation with lower order term (1.6). We supply in this paper the proof of 



elliptic concentration comparison in the two forms that are needed to try pass to the 
parabolic problem via discretization in time. However, the results are complete only in 
the case where i? is a convex function and Q = M^. Though the elliptic results are used 
here as a step towards the parabolic theory, they have an interest in themselves as an 
improvement on the symmetrization result developed in |21| . 
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• Complementing this analysis, we prove in Section [4] that one the elliptic comparison 
results that is needed to build a good parabolic theory is false in the case of a concave B 
of the form B(v) = f™, < m < 1. 

• The main issue of symmetrization for linear or nonlinear fractional parabolic equations is 
addressed in Section [5} After the iteration steps described above, the elliptic results allow 
to conclude similar comparison results for the mild solutions of the evolution problem 



(1.5) when A is concave, i.e., in the range of exponents < m < 1 of the Fractional 
Fast Diffusion Equation (FFDE), and also in the most popular case, the Fractional Heat 
Equation (FHE) 

(1.7) dtu + {-Ar/^u = f, u{0)=uo. 

A number of consequences are derived from the symmetrization result in the form of a 
priori estimates, much in the manner these consequences are derived in the case of equation 
involving the standard Laplacian. In particular, we can use the Barenblatt solutions of 
the FHE and FFDE constructed in |51j as a worst-case to obtain a priori estimates for the 
solutions in the spaces, a result that was one of the main corollaries of paper [H]. Such 
consequences are important in developing the general theory of such equations, which is 
one of the aims of the symmetrization techniques. In order to keep a reasonable length 
for this paper, we have decided to explain such consequences in a companion paper, [53] . 

• Returning to the presentation of the main results, an important gap was therefore left 
in the analysis, namely to examine what happens with this approach when applied to the 



FPME with m > 1, or more generally to (1.5) when A is not concave. To our surprise. 



the result of comparison of concentrations is false for the evolution problem, i. e., for the 
FPME with m > 1. As a consequence, we cannot use the Barenblatt solutions of the 
FPME as a worst case to obtain a priori estimates for the solutions in the spaces. The 
surprising negative results about mass comparison for parabolic equations are described 
in Section [6} This is to be seen in parallel to the counterexample found for the elliptic 



equation (1.7); by the way, this one was found later and is less intuitive. 



We conclude by a short section containing comments, extensions and open problems. 



2 Preliminaries on symmetrization 

A measurable real function / defined on is called radially symmetric {radial, for short) 
if there is a function / : [0, oo) — )• M such that f{x) = f{\x\) for all x G M^. We will 
often write /(x) = f{r), r = \x\ > for such functions by abuse of notation. We say 
that / is rearranged if it is radial, nonnegative and / is a right-continuous, non-increasing 
function of r > 0. A similar definition can be applied for real functions defined on a ball 
Br{Q) = {x e : |x| < R]. 

Now, let Q. be an open set of and / be a real measurable function on Q.. We will 
denote by |-| the A^'-dimensional Lebesgue measure. We define the distribution function 



5 



/i/ of / as 

f^f{k) = \{xGn:\f{x)\>k}\ ,k>0, 
and the decreasing rearrangement of / as 

/* (s) = sup{A; > : /i/ (A:) > s} , s G (0, \n\) . 

We may also think of extending /* as the zero function in [|r2|, oo) if is bounded. From 
this definition it turns out that nj* = fJ-f {i- e. , /, and /* are equi-distributed) and /* is 
exactly the generalized inverse of fif. Furthermore, if wtv is the measure of the unit ball 
in and f]* is the ball of centered at the origin having the same Lebesgue measure 
as fi, we define the function 

f#ix) = riu^Nlxf) ,x€n#, 

that will be called spherical decreasing rearrangement of /. From this definition it follows 
that / is rearranged if and only if / = f^. 

For an exhaustive treatment of rearrangements we refer to [7j, [31], or the appendix in 
|41| . Here, we just recall the conservation of the norms (coming from the definition of 
rearrangements and the classical Cavalieri principle): for all p G [l,oo] 



as well as the classical Hardy-Littlewood inequality (see [29] ) 



/ 1/(3 
Jn 



(2.1) / \f{x)g{x)\dx< I f*{s)g*{s)ds= I f*{x)g*{x)dx, 



where /, g are measurable functions on il.. 

• We will often deal with two-variable functions of the type 

(2.2) /: (2;,y) eCn^/(x,y) GR 

defined on the cylinder Cn := Q x (0,-|-oo), and measurable with respect to x. In that 
case, it will be convenient to define the so-called Steiner symmetrization of with respect 
to the variable x, namely the set := x (0, -|-oo) . Furthermore, we will denote by 
Hf {k, y) and /* (s, y) the distribution function and the decreasing rearrangements of ( 2.2 ), 
with respect to x for y fixed, and we will also define the function 

f# {x,y) = r{u;N\xf,y) 

which is called the Steiner symmetrization of /, with respect to the line x = 0. Clearly, 
is a spherically symmetric and decreasing function with respect to x, for any fixed y. 

• There are some interesting differentiation formulas which turn out to be very useful in 
our approach. Typically, they are used when one wants to get sharp estimates satisfied by 
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the rearrangement u* of a solution n to a certain parabolic problem, for in that context 
it becomes crucial to differentiate with respect to the extra variable y under the integral 
symbol, in the form 

J {u{x,y)>u*(s,y)} '^V 

For the sake of completeness, we recall here two formulas, of first and second order, avail- 
able in literature. The following proposition can be found in }35j . and is a generalization 
of a well-known result by Bandle (see 

Proposition 2.1 Suppose that f e H^{0,T; L'^{n)) for some T > 0. Then 

r eH\0,T;L\0,\n\)) 

and if \ {f{x,t) = f*{s,t)} \ = for a.e. {s,t) £ (0, x (0,T), the following differen- 
tiation formula holds: 

(2.3) / ^^x,y)dx= r^{T,y)dT. 

Jf{x,y)>f*{s,y) oy Jo oy 

Moreover, the following second order differentiation formula (which was also proved in 
[3] in a more regular framework) is due to Mercaldo and Ferone (see [26]): 

Proposition 2.2 Let f £ VF^'°° (Cq). Then for almost every y G (0, -|-c«) the following 
differentiation formula holds: 



f{x,y)>f(s,y) dy'^ ' dy"^ Jo ' Jf{x,y)=f''{s,y) \^==f\ 

I 



dl' ^ X 2 

y) l^^/l ' 7 \Jf{^<y)=f*is,y)\'^^f\ 



+ 1 / %^d^^-^(:r)l ( / ^dn""-' 



2.1 Mass concentration 



We will provide estimates of the solutions of our elliptic and parabolic problems in terms 
of their integrals. For that purpose, the following definition, taken from [H], is remarkably 
useful. 

Definition 2.1 Let f,g£ Lj^^iM.'^) be two radially symmetric functions on M^. We say 
that f is less concentrated than g, and we write f < g if for all R > we get 

/ f{x)dx < / g{x)dx. 
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The partial order relationship -< is called comparison of mass concentrations. Of course, 
this definition can be suitably adapted if /, g are radially symmetric and locally integrable 
functions on a ball Bji. Besides, if / and g are locally integrable on a general open set 17, 
we say that / is less concentrated than g and we write again f < g simply if /'^ -< 5*, 
but this extended definition has no use if g is not rearranged. 

The comparison of mass concentrations enjoys a nice equivalent formulation if / and g 
are rearranged, whose proof we refer to [29], |17j . |47j : 

Lemma 2.1 Let f,g £ L^{Q) be two rearranged functions on a ball Q = Bji{0). Then 
f ~< g if and only if for every convex nondecreasing function $ : [0, 00) — t- [0, 00) with 
$(0) = we have 

(2.4) / ^f{x))dx< [ Hg{x))dx. 

Jn Jn 

This result still holds if R = 00 and f,g £ Lj^^{M.^) with g as |x| —)• 00. 

Prom this Lemma it easily follows that if / -< 5 and /, g are rearranged, then 

(2.5) \\f\\Lp{n) <\\g\\Lp{n) Vpe[l,oo]. 



3 Elliptic Problems with lower order term 

3.1 Recall of existence, uniqueness and main properties 

As explained in the Introduction, the implicit time discretization scheme directly connects 



the analysis of the evolution equation (1.5) to solving the elliptic equation (1.6). There- 
fore, we start our analysis by the following nonlocal Dirichlet problem with homogeneous 
boundary condition: 



(3.1) 



{-A)''^^v + B{v) = fix) inn, 
V = on dQ, 



where is an open bounded set of R^, o" G (0, 2) and / is an integrable function defined 
in i7 (we will also take O = M^, and then no boundary condition is assumed, see below). 
We assume that the nonlinearity is given by a function B : M-|- — t- M_|_ which is smooth 
and monotone increasing with -B(O) = and B'[v) > 0. It is not essential to consider 
negative values for our main results, but the general theory can be done in that greater 
generality, just by assuming that B is extended to a function B : M_ — )• M_ by symmetry, 
B{—v) = —B(v). The fractional-Laplacian operator (— A)'^/^ acts on functions u in 
Q and is defined through the spectral decomposition of u, in terms of eigenvalues and 



8 



eigenfunctions of the Laplacian —A with homogeneous boundary conditions. Note that 
we have changed a bit the notation with respect to equation (1.6) in the introduction, by 
ehminating the constant h > 0, but the change is inessential for the comparison results. 

As explained in |15j and |14j , when working in a bounded domain the fractional Laplacian 
(— A)'^/^ can still be defined as a Dirichlet-to-Neumann map (in the same flavor of the 
construction in [16] for Q = M^), and this allows to connect nonlocal problems involving 
(— A)'^/^ to suitable degenerate-singular, local problems defined in one more space dimen- 
sion. In our case, a solution to problem (3.1) is defined as the trace of a properly defined 
Dirichlet-Neumann problem as follows. If w is a weak solution to the local problem 

-divx,y (y^^'^Vw) = in Cn, 

w = on SlCq,, 



(3.2) 



1 f/lll 

- lim yi--— (x,y)+ S(^i;(x,0)) = /(x) 



in il. 



where Cf^ := x (0, -l-oo) is the cylinder of basis BlCq, := dVL x [0, +co) is its lateral 
boundary, and is the constant (see |16j ) 



rfi 



r(f) 



then the trace of w over (7, Tv^{w) = w(-,0) =: v is said a solution to problem (3.1). To 
make this more precise, we introduce the concept of weak solution to problem (3.2). It is 
convenient to define the weighted energy space 

K'^{Cn) = jw; G H\Cn), = on ^lC^ : ^ y^-''\V ^,yw{x,y)\'^ dxdy < ooj , 



equipped with the norm 
(3.3) 

Then, following 



y^ \\7w{x,y)\'^ dxdy ] 



1/2 



we provide the following definition 



Definition 3.1 Let i7 be an open bounded set o/ and f G L^{Q). We say that w G 
^o^'^iCn) is a weak solution to (3.2) if TrQ{B{w)) =: B{w{x,0)) G L^(0) and 



(3.4) 



y^ "^x,yW ■Vx,y^dxdy + I B{w{x,0)) (f{x,0)dx = I f{x)(p{x,0)dx 



for all the test functions if G C^{Cq) vanishing on the lateral boundary 8^0^. 
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If is a solution to the "extended problem" (3.2), then the trace function v = Tr^tt; 



will be called a weak solution to problem (3.1). Concerning existence of solutions, their 



smoothness and contraction properties, we excerpt some known results from [22], [23], 
which can be extended for our more general nonlinearity B. 

Theorem 3.1 (see [22j) For any f S L°^(i7) there exists a unique weak solution w G 
XQ^'^iCn) to problem such that Trn{B{w)) G L°°(17). Moreover, 

(i) Regularity: we have w G C"(Cn) for every a < a if a < I (resp. w G C^'"(Cq) for 
every a<a — lifa>l). Arguing as in [15], higher regularity of w depends easily on 
higher regularity of f and B . 



(ii) contraction: if w^w are the solutions to (3.2) corresponding to data f,f, the 
following contraction property holds: 



(3.5) / [B{w{x,0))-B{w{x,0))]^dx< / [f{x)-f{x)]+dx. 

Jq Jn 

In particular, we have that w > in Cn whenever / > on fi. Furthermore, if we put 
u := B{w{-,0)), then for all p G [1, oo] we have 

\\u\\Lp(n) < ll/llLP(n)- 

(iii) For data f G L^(J7) the weak solution is obtained as the limit of the solutions of 
approximate problems with fn G L^{Q) H L°°(r2), fn^finL^, since then the sequence 
{B{wn{x,0))}n also converges in to some B{w{x,0), and \\B{w{x, 0))\\i < ||/||i, hence 
Vn in uniformly bounded in U' for all small p. Property (ii) holds for such limit solutions. 

We give a short account of the proof of these results for the reader's convenience. See 
more details on this issue in the forthcoming work [25]. In order to get the existence of a 
weak solution, we first define the integral function of B 

t 



G{t)= [ B{C)dC, 
Jo 



then we minimize the functional 



J{w) = — / y^-" \Vw\^dxdy + / G{\w{x,0)\)dx - / f{x)w{x,0)dx 
2Kcr Jcn Jn Jn 

over the space Xq^'^{Cq). where A := B^^ is the inverse of B. Arguing as in [23], and 
using the Holder, trace, and Young inequalities we find that the functional J is coercive 
on X. In order to prove that J' is weak lower semi-continuous, let be a sequence 

in Xq^^(Cq) converging weakly to vu. By the trace embedding theorem, we have (up to 
subsequences) 

Wni-, 0) w{-,0) Strong in L'^{n) Wq G [1, 2N/{N - a)), 
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then 



Wn{-,0) ^ w{-,0) a.e.inVL. 



Now Fatou's lemma implies that J is weakly lower semicontinuous on Xq^^(Cq). Then 
there exists a minimizer w £ Xq^'^{Cq) of J^. Furthermore, a truncation argument shows 
that we can suppose 'w{-,0) G L°°(il) and 

lk(-,o)||i^(n) <^(||/IU^(Q)). 

Finally, computing the first variation of in the direction of any </? E XQ^'^{Cfi) we obtain 
that w is a weak solution to (3.2) in the sense of definition (3.1 ). The contraction property 



(3.5) in Theorem 3.1 follows by the arguments of |22j-|23j. 

• Let us now consider problem (3.1) in the whole M^, where the fractional Laplacian is 



defined by a singular integral. The problem is 

{-A^^^v + Biv) = fix) 

v{x) — )• 



(3.6) 



in M- 



as X 



N 



OO, 



where / E L\R^) D L°°(M^), and we can define again a suitable meaning of weak 

solution, making use of a proper extension problem. Indeed, if we denote by (Cj^jv ) , 
the upper half-space, the completion of C°°(Cjj 

pN 



being 



with respect to the 

norm (3.3) with Q, replaced by M^^, then a solution v to (3.6) is the trace on of a weak 

to the problem 



solution w E X''/'^{C 
div 



(3.7) 



x,y 



'Vw) 







m 



X (0,+oo) 



1 Quj 
— Vim y^-" —{x,y)+ B{w{x,G)) 
Ka y^o+ dy 



fix) 



X E 



Of course, we mean that w E X"^'^iC-g_N) is a weak (energy) solution to (3.7) if equality 



(3.4) holds, with replaced by M . In order to obtain the existence and uniqueness of 
solution to problem (3.7), we can ague as in |23]. For any i? > 0, we consider the solution 
wr to (3.2) corresponding to the data fR = fxBuio)! where Q is the ball i?_R(0) centered 



at the origin. If the data / is nonnegative (the case of changing sign data can be treated 
as in [22]), we obtain an increasing sequence of non- negative solutions {wr} converging 



to a weak solution w to the problem (3.7) in the upper half-space. Then the contraction 
property (3.5) holds in M^, from which uniqueness and preserving sign property follow. 
i?(u;(-,0)) then we have 



Moreover, if u 



< 



\u\ 



< 



Remarks. The approximation method we have used to prove the comparison theorem 
in the whole actually says that we can approximate the solution v to problem (3.6) 
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with the fractional Laplacian on M by a sequence of solutions of Dirichlet problems of 
the type (3.1) with the fractional laplacian defined on balls, with homogeneous boundary 
data. 



We point out that Theorem |3.1| and the related considerations of existence of solutions 
on R-^ still hold if B : M — )• M is assumed to be increasing and B{0) > (and this remark 
will enter in Subsection 3.5). If we want to extend B to the whole real axis it suffices to 
set B{-v) = 2B{0) - B{v) for all v > 0. 

From now on, we will always assume that the right-hand side / is nonnegative. 



3.2 The extended problem and concentration comparison 



Let us address the comparison issue. Our goal here is to compare the solution v to (3.1) 
with the solution V to the problem 



(3.8) 



{-/^yl^V + B{V) = f*{x) in 9.* 

y = on an*. 



A reasonable way to do that is to compare the solution w to (3.2) with the solution to 
the problem 



(3.9) 



' - dW^,y (yi-'^VV') = 
= 



in Cq# 
on OlC^^ 



1 



lim y^-'' ^(x, y) + B(Mx, 0)) = f* (x) in n*, 

y^o+ ay 



where V'(3;,0) = V{x). According to [21j, using the change of variables 



problems (3.2) and (3.9) become respectively 



(3.10) 



— z 

dw 



A^w = 



dz 



(x,0) 



a- 



cr-l 



if {x) - B{w{x,0))) 



in Cq 
on dLCa 

in 



12 



and 



(3.11) 



9> . 



dip 



■ n# 



dz 



on OlC* 

{x, 0) = a^-^K^ (/# (x) - B{iP{x, 0))) m 0#. 



where ly := 2(cr — l)/(7. Then, the problem reduces to prove the concentration com- 
parison between the solutions w{x,z) and ip{x,z) to (3.10)-(3.11 ) respectively. Fohowing 
[21j, using standard symmetrization tools (among which the differentiation formulas in 
Propositions (2.1)-(2.2) are essential), if we introduce the function 

(3.12) Z{s,z)= r{w*{T,z)-r{r,z))dT, 

Jo 

then we get the inequality 

(3.13) -z''Z,,-p{s)Zss<0 

for a.e. (s, z) £ D := (0, |il|) x (0, +00) . Obviously, we have 

(3.14) Z{0,y) = Z,{\n\,y)=0. 

A crucial point in our arguments below is played by the derivative of Z with respect to 
z. Due to the boundary conditions contained in (3.10)-(3.11 ), we have 

(3.15) Z,is,O)>0^ [\b{w*{t,0))- B{r{T,0))dT 



where 

Now observe that the function 

Y{s,0) 



B{w*{t,0)) - B{llj*{T,0))dT 



has the same points of maximum or minimum and the same regions of monotonicity than 
Z{s,0). 



3.3 Comparison result for concave B 



Theorem 3.2 Let v be the nonnegative solution of problem (3.1) posed in a bounded do- 
main with zero Dirichlet boundary condition, nonnegative data f E L^{U) and nonlinearity 
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B(v) given by a concave function with B(0) = and B'{v) > for all v > 0. IfVis the 
solution of the corresponding symmetrized problem, we have 



(3.16) 

The same is true if 



t;*(x) -< V{x). 



oN 



Proof. In this case we pose the problem first in a bounded domain Q of with smooth 
boundary. We also assume that / is smooth, bounded and compactly supported, since 
the comparison result for general data can be obtained later by approximation using the 
dependence of the map / i— t- B{v). 

(i) We want to prove that Z(s,0) < for all s G [0, It is easy to prove that a 
positive maximum of Z{s, z) cannot happen at the lateral boundary s = |0| for z > by 
the stated boundary conditions and the Hopf's boundary principle. 

In order to study the possible positive maximum at the line z = we proceed as follows. 
The concavity of B implies that for a, & > we have B{a) — B(b) > B'{a){a — b). Using 



this and (3.15), it follows that 



Z,is,0) > 0.y(s,O) > 9^ / B'iv*{T))[w*{T,0)-rir,0)]dT 

Jo 

B'{v*{T))Zs{T,0)dT. 



(3.17) 



If we set g{s) := B'(v*{s)), since B' is decreasing, we notice that g is an increasing function 
bounded from below by g{0) = B'{\\v\\oo)- If the positive maximum of Z{s, 0) happens for 



s = So then, using an integration by parts in (3.17) we can write 

rsQ 



(3.18) 



Z^{so,0)>9a B'{v*{s))Zs{s,0)ds 

'0 

rso 



9n 



5(0)Z(so,0) + r\ziso,0) - Zis,0)]dg{s) 

JO 



> 



which is impossible because ^^(50,0) < 0. Then Z(so,0) < 0, that is Z < 0, namely 



W*{T,z)dT< / ip*{T,z)dT. 

Jo 



Another remark is that either Z = or 



(3.19) 



Z < in (0, |17|) X [0,oo) : 



indeed, if Z ^ for the previous arguments it cannot be Z = in some points of 
(0, X (0,oo) (otherwise it would reach the maximum in this domain, hence it would 
be constantly by the maximum principle). On the other hand, if Z{so, 0) = for some 
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point So £ (0, by the Hopf boundary maximum principle we have Zz{sq, 0) < 0, but 



by (3.18) we have Zz{so,0) > 0. 



(ii) Here is a simpler proof in the important special case of the linear fractional diffusion, 



i.e., when B{v) = v. Indeed, from (3.15) we have the inequality 

z,{s,o) > e,z{s,o). 



Now we simply observe that (3.13) can be rewritten as 

-p{s)-'Z^ 



ZZ Z Zgg ^ 



therefore, multiplying both sides by Z+ and integrating by parts over the strip [0, |f2|] X 



(0, +oo), the boundary conditions (3.14) and the fact that Z{s, z) — )• as z — )• oo imply 



oo p\n\ 



piui poo p 

/ p{s)-'Z,{s,0)Z+{s,0)ds+ / / 
Jo Jo JO 

+ / / p{s)-^\{Z+)J^dsdz <0 
Jo JO 



z-n (z. 



ds dz 



namely 



/•oo j-oo 

/ / z-''\{Z+)^\^dsdz+ / / p(s)-^\{Z+)^\^dsdz<Q 
Jo Jo Jo Jo 



oo p\n\ 



hence = 0. 

• Problem in the whole space. The previous arguments still apply if the problem is posed 



in the whole space M , namely if v solves (3.6). Indeed, in this case we may use the 



boundary condition Zs{s,y) — )• as s — )• oo. Alternatively, according to what remarked 
in Section 3, we may approximate the solution w to the elliptic problem in the upper 
half-space (3.7), with nonnegative / G L^(M^)nL°°(M^), with the family wr of solutions 



to problems of the type 

' - div^,y {v^-^Vwr) = 

WR = 



(3.20) 



in Cbh, 
on OlCbh, 



1 dzL) 

lira y^''' ——{x,y) + B{wr{x,0)) = fR{x) in Br, 

Ka y'^o+ oy 



where Br is a ball of radius R at the origin. According to Theorem 3.2, we obtain 



(3.21) 



w*R{T,y)dT < / i;*ji{T,y)dT 
Jo 
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for all s G [0, and y > where ipji is the solution to 
' - div^,^ {y^-'^VijR) = 



in Cbr 
on dLCsn 



-— hm i?(^/>,j(x,0)) = /#(x) in Br. 

Then we get (see Theorem 7.3 in [23]) wr — )• w and iI)r — )• V' as i? — )• oo, where iIj solves 
- div^,j^ (yi-'^V^) =0 m X (0, +oo) 



(3.22) 



1 dib 
lim y^-'' -^(x, y) + Bi^ix, 0)) = f* ix) x G M^. 



Therefore, letting i? — )• oo in (3.21) we find 



w*{T,y)dT < / i;*{T,y)dT 

JO 



for all s > and y > 0. 

Remark. We also wanted to prove that Y{s, 0) < 0, i. e., 

(3.23) / B{v*{T))dT< / B{V*{T))dT. 

Jo Jo 

but it did not work. See next section. 



3.4 Comparison of concentrations for radial problems 

This second result is a variation and extension of the previous comparison result. We 
consider the same assumptions on B and O. 



Theorem 3.3 Let vi,V2 be two nonnegative solutions of problem (3.1) posed in a ball 



Br{0), with R G (0, +oo] with zero Dirichlet boundary conditions if R < +oo, nonnegative 
radially symmetric decreasing data /i,/2 G L^{Br{0)) and nonlinearity B{v) given by a 
concave function for v >Q, with B(0) = and B'{v) > for all v > 0. Then vi and V2 
are rearranged, and 



(3.24) 



/i ^ /2 implies vi < V2 ^ 



Proof. As in the proof of Theorem 3.2 we arrive at the inequality 
(3.25) -z-Zl'^-p{s)Zlf<Q 
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satisfied a.e. in the strip (0, |i?ij(0)|) x (0, +00) by the function 

Z^'^{s,z) = [ {wl-W2)dT 



where wi and W2 are the solutions of the extensions problems associated to vi and V2 
respectively. Concerning the boundary conditions, since we have /i -< f2 we get 

Zl'\s,0)>9, riB{wl{T,0))-B{w*2{T,0))dT+ fir^-fDdr 
Jo Jo 

> [ {B{wI{t,0)) - B{w2{T,0))dT. 







Then we conclude as in the proof of Theorem 3.2 



3.5 Comparison results for convex B 

In the case of a convex nonlinearity we prove a stronger result in the whole space. For the 
sake of clarity, we first prove the result when B is a superlinear nonlinearity in the sense 
that is made precise next: 



Theorem 3.4 Let v be the nonnegative solution of problem (3.6) posed in Q = M^, 
nonnegative data f G L^(M^) and nonlinearity given by a convex function B : M+ — ?■ M+ 
which is smooth, and superlinear: B{v) > ev for some e > and all v > 0. If V is the 
solution of the corresponding symmetrized problem, we have 

(3.26) V* -< V, B{v*) -< B{V). 

Remark. The simplest example of superlinear nonlinearity is of course the linear case, 
B{v) = cv. A nontrivial example from the literature would be the remarkable nonlinearity 
A{t) = log(l+t) in the model of logarithmic diffusion, ^^4^. Then, B(s) = A~^{s) = e* — 1, 
s > 0. We will relax the restriction B{v) > ev below by approximation. 

Proof. In order to prove that Z(s, 0) < for all s G [0, 00) we argue as follows. We have 
B{w*{t,0)) - B{ip*{T,0) = B'{^){w*{t,0)) - V*(-r,0)), where C is an intermediate value 
between w*{t,0) and ^l^*{T,0). Since B is convex, B' is an increasing real function and 

B{w*iT, 0)) - Birir, 0)) < B'{w*iT, 0)))(ii;*(r, 0)) - (r, 0)) 



Due to the maximum principle and the boundary conditions (3.14), unless Z is constant, 
the maximum of Z can be achieved either on the half-line {(0, z) : z > 0} or on the segment 
line {(s,0) : s G [0,oo)}. Suppose this second circumstance occurs, and let (so,0) be a 
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maximum point. Assume sq > 0. We also have Zz{sq,0) < by Hopf's maximum 
principle, and by (3.15), this leads to Y{so,0) < 0. Then for s > sq 



Y{s,0)-Y{so,0)= / [B{v*{T))-B{V*{T))]dT 

s " 
< / B'{w*{T,0))iv*iT)-V*iT))dT. 



so 



After integration by parts 

Yis,0)-Y{so,0) < [B'{v*iT))iZiT,0)-Z{so,0))]: 
B"iv*iT))v:iT){Z{T,0) - Z{so,0))dT. 



Since Z has a maximum at sq and B' is positive, the first term in the RHS is nonpositive. 
As for the second, we have: B" > 0, v* < 0, and Z{s,0) — Z{so,0) < 0, hence the last 
term is also nonpositive. We conclude that Y(s,0) < Y{so,0) < for all s > sq. This is 



a contradiction, because by the conservation of mass property (see proposition 3.5 below) 
we have Y{oo, 0) = 0. Then sq = and Z < 0. 

(ii) Once we have Z{s, z) < we also want to prove that Y{s, 0) < 0. We use the fact 
that s = is a point of maximum of Z and write 

Yis,0)<[B\v*{r))ZiT,0)]l- f B" {v* {t))vI{t)Z{t,^)) dr < ^. 





Also, we obtain the same result by using Lemma [2.1[ taking advantage of the convexity of 
B and choosing any convex, increasing function $ : [0, oo) — )• [0, oo). This ends the proof 
of the concentration comparison theorem in this case. □ 

• The only remaining question here is then to prove that ||i?(u;(x, 0))||ii = ||-B('(/'(a;, 0))||ii. 
Under the additional assumption B{s) > es for all s > 0, this will be a consequence of 
the following mass conservation result for the solutions of the elliptic equation. 

Proposition 3.5 Let v be the weak solution of {-A^/'^v + B{v) = f with f E L^(R^) 
nonnegative and let u = B{v), with B satisfying the same assumptions as in Theorem 3.4 
Then we have 

u{x) dx = f{x) dx. 



Proof. Using a nonnegative nonincreasing cutoff function C(s) such that C{s) = 1 for 
< s < 1 and = for s > 2, we rescale such function to Cr{x) = ({\x\/R). Then we 
have 

(3.27) / f{x)CR(.x)dx- [ u{x)CR{x)dx= [ v{{-Ar/^CR)dx 

JR^ JR^ JR^ 
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Due to the superlinearity assumption, if u{x) = B{v{x)) we get 



1 



c 



< - / \u{x){-Ar/^CR\dx< 



\u{x) \ dx 



which in the hmit — )• oo tends to zero. □ 

• The property of mass conservation for solutions in the whole space is also true for 
some convex B that are not superlinear, like B{v) = with some p > 1 but near 1, 
but it is not true for all p > 1. However, the comparison result we are looking for 



(which will extend Theorem 3.4 for such kind of nonlinearity) will be true and the proof 
proceeds by an approximation process, approximating B{v) by B^iv) = B{v) + ev, solving 
the approximate problem, deriving the comparison result and passing to the limit. The 
details are as follows. 



is smooth, convex, B(0) = 



Proposition 3.6 Suppose f E L^(R^) and B : ^ 
and B'{v) > for all v > 0. Let Ve he the solution of problem (3.6), with nonlinearity 
given by Be{v) = B{v) + ev. Then ^ v as s ^ pointwise and in L^(M^). 



Proof. We first prove the result on a bounded domain Cl. Suppose that / G let 
Vs be the solution to (3.1), and let be its a— harmonic extension to the cylinder Cn. 
We have that Wr is obtained as minimizer of the functional 



Je{w) 



2k„ 



y^-''\Vw\'^dxdy+ I Ge{\w{x,0)\)dx - I fw{x,0)dx 



with 



Ge{t) 



[B{C) + eC] di 



over the space Xq^'^{Cq). Moreover the trace Ve over Q. of We is bounded and 
(3.28) \\we{-M\L^m < AiWfUoo^^^. 



Taking Ws as a test function in the weak formulation of problem (3.2), namely in the 
formula 
(3.29) 



,-"V...V^<!x.,+ /BK(x,0))„(.,0)d.+. /'..(..0)^(x,0)d. = ., / 

Jo. Jn Jn 



f{x) (p{x. 



for if € Xq^'^{Cq), the Young and trace inequalities imply that {ws} is bounded in 

Xq^^(Cq). Then we can extract a subsequence {we} (we used the same labeling for 
simplicity) such that 

We ^ w weak in Xq^'^{Cq). 
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Then the compactness of the trace embedding inequaUty gives 

We{-,0) ^ w{-,0) ase^O strong in L''(J7) £ [1,2N/{N - a)). 



Using (|3.28p, Lebesgue's dominated convergence imphes 

B{we{-,0))ip{-,0)dx / B{w{-,0))ip{-,0)dx 



for all if G Xq^^(Cq). This is enough to pass to the limit in (3.29) and obtain that w is 



the weak solution to (3.2), that is w = w{-,0) solves (3.1) 

If / e L°°(R^), let V be the solution to ( [3^ . We know that the solution to ( [3^ ) 
with nonlinearity is the trace on of the solution We to the problem (3.7), with 



the nonlinearity B^. By the arguments we recalled in Subsection 3.1, we have that the 
sequence of solutions {w^^j^^g to problem (3.2), defined on the cylinder C^^(q-^, with 
nonlinearity B^ and data = /xb^(o)) converges pointwise to as ii oo, that is 



(3.30) 



R 



pointwise as i? — )■ oo. 



Moreover, if = w^{-, 0), the arguments explained above show that 



(3.31) ^ as e ^ strong in L'^{Br{0)) \/q € [1, 2iV/(iV - a)) 

where is the solution to the Dirichlet problem (3.1), posed on the ball Bji{0). In 
addition, we have 

(3.32) f — )• v pointwise as ii —)• oo. 

Then (3. 30), (3. 31), (3. 32) implies that converges to v pointwise and in L-'^(]R^). □ 



Then we are able to prove the comparison result of Theorem 3.4 also for power nonlin- 
earities like B{v) = for large p > 1 (which means that we can include the fast diffusion 
range when we pass to the parabolic setting in Section [5]) 



Theorem 3.7 Let v be the nonnegative solution of problem (3.6), nonnegative data f G 
L^(M.^) and the nonlinearity given by a convex function B : M+ — t- M+, with B{0) = and 
B'{v) > for all V > 0. If V is the solution of the corresponding symmetrized problem, 



the conclusion of Theorem 3.4 still holds. 



Proof. By virtue of Theorem 3.4 we have 

(3.33) V* -< V„ Bivf) -< BiVe) 

for all e > 0, being v^, the solution of problem (|3.6|) and its symmetrized with the 



nonlinearity Bf,. By Proposition 3.6 we have that for all s > 



vZ dr 



V* dr, 



v:dT 



V* dr 



as e — )■ 0. Passing to the limit in (3.21) we find the desired result. 



□ 
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3.6 Second comparison result for convex B 



Here is the second result, about comparison of concentrations for radial problems. We 
leave the proof to the reader. 



Theorem 3.8 Let vi,V2 be two nonnegative solutions of problem (3.1) posed mQ. = R^, 
with nonnegative radially symmetric decreasing data /i, /2 G L^{^) and nonlinearity B{v) 
given by a convex function with B{0) = and B'{v) > for all v > 0. Then, vi and V2 
are rearranged, and for fi -< /2 we have 

(3.34) vi{x) ^ V2{x), B{vi) -< B{v2) ■ 

Remark. These results are in perfect agreement with the results of |45J for the standard 
Laplacian case. 



4 Counterexample for elliptic concentration comparison with 
convex powers 

If we compare the results of the preceding section for convex B and concave B we realize 
that the conclusion is weaker in the latter case. This seemed to us a possible defect in the 
technique since in the case of the standard diffusion a = 2 the results are identical. But it 
turned out that in the fractional equation the concave case has an essential difficulty. 

Our aim is now to prove that actually the general concentration comparison does not hold 
in the whole M^, for an equation of the form 

(4.1) hL^{u"') + u = f, L^ = {-Ar/\h>0, 



with m > 1. The precise result is stated in Theorem 4.3 The reduction to power-like 



nonlinearity simplifies the calculations and is the most important case in the applications. 



The solutions satisfy u — )• as |x| — )• oo. Comparing this equation to (3.1), we notice 
that here we denote hy u = B(v) = v^l'^ the unknown function. Equation (4.1) is posed 
in with nonnegative and integrable data /(x). We want to describe the asymptotic 
behaviour of the solution as |x| — ?■ oo, more precisely its rate of decay. We will focus on 
the dependence of the behaviour on the constant /i > 0. This constant is important since 
it represents the time increment when discretizing the evolution problem. We stress the 
dependence by often denoting the solution as u{x;h). Our main result says that roughly 
speaking 

(4.2) u{x;h) ~ h\x\-^^+^'^ 
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when is large and h small. We assume m > 1. Note that the parameter h can be 
changed, or fixed to 1, by using the scaling 

(4.3) u{x) = au{bx), f{x) = af{bx). 

If a^~^b'^h = h, then u satisfies: /iLo- (5™) + u = f. In other words, u{x; h) = au{bx; h). 
This will be of great use in deriving the negative implication for the concentration analysis 
in Section l43l 



4.1 Elliptic positivity estimate via subsolutions 



The first step in our asymptotic positivity analysis of solutions of (4.1) is to ensure that 
solutions with positive data remain positive in some region. We only need a special case 
that we establish next. 



Lemma 4.1 Let u{x;h) be the solution of (4.1) with RHS f{x) > such that f{x) > 1 



for \x\ < 1. We assume that m > 0. Then, for every R < 1 there are constants Ai,hi > 
(depending on R) such that 

(4.4) u{x;h)>Ai for \x\ < R, < h < hi. 

Proof. • First, we construct a subsolution for a related problem that has an explicit form 
and compact support. Let g{x) be the explicit function, 

g{x) = ^{1 - r%/\ r=|rE|>0. 

Getoor [2H], Theorem 5.2, proves that L(jg[x) = cq > on the ball of radius 1 where 
g is positive, while ^0-/25 < for r > 1, with an explicit formula that goes to minus 
infinity as r — t- 1 and behaves as ~ j.-{n+(t) -^j^g,^ go. Next, we consider the following 
combination 

/i(x) ■.= hL„g + g^/^. 



This can be seen as follows: the solution u of equation (4.1) corresponding to RHS /i is 
ui = g^/"^. 

Let us now try to estimate /i: for r > 1 we have /i < 0. For r < 1 we have /i = 
Coh + > 0, besides /i < (1/2)^"^ + /icq < 1 if /i < (1 - 2-")/co. Under these 

restrictions on h, ui = g^^'^, the solution for RHS /i, serves as a subsolution for the 
RHS f{x) = Xi(0)) the characteristic function of the ball of radius 1. This means that 
the solution u corresponding to such / is equal or larger than ui = g^^^. Since ui{x) is 
uniformly positive in the ball of radius 1/2, u{x) is uniformly positive in the ball of radius 
1 /2 when < h < hi. 

• By means of scalings to put the dimensions in x, u and h as in the statement, n 
We now proceed with the asymptotic estimate from below. 
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Theorem 4.1 Let u{x; h) he the solution of (4.1) with RHS f{x) > such that f{x) > 1 



in the ball -Bi(O). We assume that m > 1. Then there are constants C-,Ri,hi > such 
that 

(4.5) u{x;h)>C-h\x\-^^+''^ 
if \x\ > Ri and < h < hi. 



Proof. We will use the standard comparison theorem to reduce the case where / is a 
smooth version of the characteristic function of the ball -Bi(O). Then known theory says 
that u < 1 everywhere and is continuous, radially symmetric and decreasing in r = \x\. 
In fact, a bootstrap argument shows that u G C°°. Since ||u||i < ||/||i = con we also have 
a first decay for u near infinity of the form 

n(r) < Cr-^. 

Of course, this first estimate is not sharp, in view of our next results. 
• We want to construct a subsolution of the form 

(4.6) [/™(x;/i) = G(|x|) + /i™F™(|x|), 

which will be valid for < h < hi. Here the functions G,F > and the constant hi > 
have to chosen carefully, as explained below. 

We take G{r) = for r = |x| > 1/2 so that U{x; h) = hF{x) there. If G is also smooth 
we have L^j G bounded and we can also choose G so that 

LaG < -Cir-(^+<^) for r > 1/2. 

We may choose as G a smoothed version of the previous Getoor function, using convolu- 
tion. 

We also need F positive, smooth and F{r) ~ C2?"^(^+'^) as r — )■ oo to get the desired 
conclusion after the comparison argument: u{x; h) > U{x; h) > C hr~^'^~^'^^ (if r is large 
and /i ~ 0, see below) 

To check the property of subsolution we proceed as follows. We have 

L^U"" = L„G{x) + h"'LaF'''{x) 

As we have pointed out, our choice of G leads to the above estimate for L„G with negative 
sign. We also have F < G2r~^'^~^°'^ for r > 1/2, by Lemma 2.1 in [T3] we have that since 
^ (9(^-{iV+a)m) ^^^^ _^ > ^ ^^^^ choosc F SO that IL^F"*] < C3r-(^+'^) 

for some positive constant C3 and r > 1/2. Then we will have 

U + hL^U'^ <h{F + L^G + h^'L^F"') < h{C2r-^^+^^ - Cir-^^+^'> + h'^'L^F""), 
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which will be negative for all r > 1/2 if 

C2 + hJ^C-i < Ci . 

In order to make sure that such constants can be obtained, we fix first G and this de- 
termines Ci. We then use a tentative Fq for the function F and multiply it by a small 
constant so that F and smaller than C\j1. Indeed, we can take C2 < Ci/2 

and h < (Ci/2C3)i/™ =: hi. Finally, we may take hi = 1. 

• The next step is to use the viscosity method to compare u and U in the Q = {\x\ > 1/2}, 
and this will prove that U{x; h) < u{x; h) in Q if h < hi. 

The following inequality establishes a suitable comparison of the boundary conditions at 
|x| = 1/2: 

U{x; h) = /iF(l/2) < ^1 < u{x; h). 



Here we use Lemma (4.1) with the choice R = 1/2, which gives u{x;h) > Ai for some 



constant Ai, \x\ < 1/2 and h sufficiently small. Now 

hF{l/2) < 2^+^ /11C2 < Ai 

up to choose C2 properly and h under a further bound. Once this is justified, we argue at 
the first point where u and U touch. Actually, we must use and approximation instead 
of u. □ 

Remark. The only restriction on m is m{N+a) > N, which means m > mi = N/{N+cr). 
4.2 Upper bound estimate 



Theorem 4.2 Let u{x;h) be the solution of (4.1) with RHS f such that < f{x) < 1 in 
the ball -Bi(O) and f{x) = for \x\ > 1. We assume that m > 1. Then there are constants 
C+,R2,h2 > such that 

(4.7) u{x;h) <C+h\x\-^^+''^ 
if \x\ > i?2 and < h < h2. 

Proof. We will construct a super-solution of the form 

(4.8) U"'{x;h) = G{x) + b"'h'^F{x)'^, 

where F > is chosen as before and 6 > 0. We will use the fact that F{r) ~ Cir^^^'^'^^ 
as r —7- 00. It follows that there is a large constant A; > such that 

kF + L^F"" >0 everywhere in . 
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Next we choose G > compactly supported in a ball of radius i?i > 1, and such that 
LfjG = Co > on the support. As r — )• oo, we get the usual L^jG ~ —C r~(^+°"). We also 
need G{1) > 1. Note that for G = we have U = bhF. In any case U is nonnegative, 
f7 > 0. We also have 

L^U'^ = L^G + hVL^F'^ix) 
We perform an analysis by regions. Thus, when G = we have 

U + /iL^f/™ = h{hF + V^K'^L^F'^ + L^G) > 

The final inequality is obtained as follows: we first put b > bo so that 

ib/2)F + L„G>0 

(recah that L^G = 0{r'^^ +"'>)). Then we put b > 2k{bh)"' to have 

{b/2)F + {bhri^F"" > 0, 

i.e. if b"^~^h'^k < 1/2; this imposes an upper bound on h. 
On the other hand, where G > we have 

U + hL^U"" > h{cQ + {bhy^L^F"") > 

if G^{bh)'^ < Co (we use the fact that L^F^ is bounded). Both conditions are fulfilled if 
< /i < /l2. 

• Now the viscosity method works in the region Q = {\x\ > 1}, and this will prove that 
U{x,h) > u{x,h) in Q. Indeed, the boundary condition at r = 1 is 

U{1) > G(l) > 1 > u{x;h). 

by the maximum principle. This ends the proof. □ 

Remark. The only restriction on m is m{N+a) > N, which means m > rrii = N/{N+a). 



4.3 Scaled data. Negative concentration result 

We have done the argument for a solution with data of height 1 supported in the ball of 
radius R = 1. If we want to change the radius to i? 7^ 1 and the height to A we can use 
the scaling 

(4.9) u{x; h) = Au{x/R; A'^-^R-'^h) 

Using this formula and the result of Theorem |4 . 1 1 applied to u, we see that the comparison 
result holds in an /i-interval of the form 

< /i< hi{R) = 
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and the new result is 



(4.10) 



u{x; h) > C- 



valid for large x and h suitably small. We are interested in conservation of mass, i.e., 
A = R-^. In that case, denoting the new solution by uji{x; h) we have 



(4.11) 



ur{x; h) > 



C- 



h 



Theorem 
(4.12) 



4.2 



In the same way, the scaling formula (4.9) applies in combination with the result of 
in the /i-interval: < /i < h2{R) = /i2 RA~^"^~^\ and the new result is 



u{x; h) < C+ 



hA'^R 

N+o 



\X 



Under conservation of mass, A = R , denoting the solution by uii{x] h) we have 



(4.13) 



ur{x; h) < 



h 



JlN{m-l) \x\N+^' 

We are ready to arrive at a contradiction in the comparison of concentrations. 

Theorem 4.3 If m > 1 there exist two nonnegative, compactly supported, bounded, radi- 
ally symmetric and rearranged functions, f and f^, such that fn -< f, and nevertheless 
the corresponding solutions u{x; h) and ur{x] h) do not obey the same relation. 



Proof. Let us choose / = xii then rescaled function fR{x) = R~ f{x/R) = R^ XR is 
compactly supported in the ball Br, has height R~^ and it less concentrated than / if 
R> 1. Let us consider the solutions u and ur that they produce, with the same coefficient 



h. If we apply Theorem 4.1 and inequality (4.13) to u and ur respectively, we have that 
ur{x) < u{x) if C-f < C-R^^'^~^\ on the condition that x is large enough, and h is small 
enough: 

\x\ > R-, \x\ > R+R; h<hi, h< /i2i?^('"-^)+'". 



Since ur and u have the same mass (because they are solutions to equation (4.1) with 
the same h and data having the same mass, see the remark below), this means that ur 
cannot be less concentrated than u for such small values of h. 

Once we have the contradiction for the equation with some h we may put /i = 1 by 
scaling. □ 

Remark. Here we see that the contradiction is obtained only for m > 1. For m < 1, C+ 
will always be larger than C-R^^'^~^\ and there is contradiction, just as predicted by 



the theory, cf. Theorem 3.8 
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Remark. In order to prove the conservation of the mass for the nonhnearity A{u) = 
with m > 1, we can argue as in Proposition |3.5[ Indeed, suppose that v be the weak 
solution of (— A)*^/^!; + B{v) = f with / G L^(R^) nonnegative and let u = B{v), with 
B : M_|_ —7- M-|_ be a concave function, strictly increasing, such that -B(O) = 0. Suppose 
first that / G L°°(M^) and |/| < K. By Theorem 
the convexity oi A = B^^, we have that the function t E 
then 



we have that u = B(v) < K. By 
— )• A(t)/t G M+ is increasing, 



V 

u 



A{B{v)) ^ A{K) 



B{v) 



K 



Now if C is the usual cutoff function and C,r{.x) = Ci^x) is its rescaled version, we still 



find equation (3.27). We also have 

A{K) 



< 



K 



\u{x){-/\r'\R\dx< 



\u{x)\ dx 



which in the limit — )• oo tends to zero. Then by (3.27) we conclude 



u(x) dx 



f{x) dx. 



If / is in L^(K^) we proceed by approximation. 

Last Remark. We want to point out the comparison performed in this section looks too 
contrived. Actually, this is partly due to the fact that it is the translation into the elliptic 
framework of the more natural parabolic counterexample, constructed in Section [6j 



5 Symmetrization for the parabolic problem 



For simplicity of exposition, we start with the case / = 0. We briefly remind the concept 
of weak solution to the Cauchy problem associated to a fractional parabolic equation : 

[ ut + {-Ay/^A{u) = x£R^,t>0 

(5-1) ^ . ^ 

I u{x, 0) = ^0(2;) X G 



Here, uq is an integrable function on M^, nonnegative in our applications), the nonhnearity 
A{u) is a nonnegative concave function with ^(0) = and A'(u) > for all u > 
(extended antisymmetrically in the general two-signed theory). Set B := A~^. As in the 



elliptic case we rewrite problem (5.1) as the following quasi-stationary problem 
- div^,s, {y^^'^Vw) = {x, y) G X (0, 00) , t > 



(5.2) 



1 



lim y 



i-o- 



dw 



= 
dB{w) 



3/-s>o+ " dy dt 
^ «;(a;,0,0) = A{uq{x)) 







X G 
X G 



pN 



y = 0,t>0. 



Then we have the following definition 
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Definition 5.1 We say that w is a weak energy solution to problem (5.2) ifw G Lf^^{(0, oo);X"'^'^ 
the function u{x,t) := {B{w{x,0,t)) is in the space C([0, oo); L^{W^ and the following 
identity holds 



^ f dip , , 
/ u—dxat 

jRiV dt K. 



'"■o- Jo JCjjjv 



for all test functions ip £ Cq(M_,_ x [0,oo)). Finally, the initial data are taken in the 
sense that 



limu(-,t) = uo(x) G LVM'M. 
t-^o 



If is a solution to 
the Cauchy problem (5.1 
uniqueness of weak so^ 



5.2), we sill say that u{x,t) : 
We refer to El 



{B{w{x,0,t)) is a weak solution to 
for questions related to existence and 



utions to problem (5.1). 



The theorems we are going to prove, Theorems 5.3 and 5.4 will come from the com- 



bination of two ingredients: the existence of a mild solution to problem (5.2) that is 



reduced to solving some elliptic problems by applying the Crandall-Liggett theory for 

n already proved for elliptic 



m-accretive operators, and the comparison theorems 3.7 



problems. Therefore, we will devote a subsection to review this material for the reader's 
convenience. 



5.1 Abstract evolution equations and accretive operators. The semi- 
group approach 

Let X be a Banach space and A : D(A) C X ^ X a nonlinear operator defined on a 
suitable subset of X. Let us consider the problem 



(5.1) 



u'{t) +A{u) = f, t>0, 
n(0) = no , 



where uq £ X and / G L^{I', X) for some interval / of the real axis. For a wide class of 
operators, in particular the ones considered in this paper, a very efficient way to approach 
such problem is to use an implicit time discretization scheme that we describe next. 
Suppose to be specific that I = [0, T] (but this can be replaced by any interval [a, b] and 
the procedure is similar). The method consists in taking first a partition of the interval, 
say, tk = kh for k = 0, 1, ... n and h = T/n, and then solving the system of difference 
relations 



(5-2) + A{uh,k) = fl 



for A; = 0, 1, . . . n, where we pose Uhfl = uq. The data set j/^^'' : A; = 1, . . . , n| is supposed 
to be a suitable discretization of the source term /, corresponding to the time discretization 
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we choose. This process is cahed implicit time discretization scheme (ITD for short) of 
the equation u'{t) + A{u) = f. It can be rephrased in the form 

Uh,k = Jh{uh,k-i + hfl^'') 

where the operator 

Jx = {I + XA)-\ A > 

is called the resolvent operator, being / the identity operator. Therefore, the application 
of the method needs the operator A to have a well-defined family of resolvents with good 
properties. When the ITD is solved, we construct a discrete approximate solution {it/i,fc}^. 
By piecing together the values u^^k we form a piecewise constant function, Ufi{t), typically 
defined through 

(5.3) Uh{t) = uh,k ift£[{k-l)h,kh] 

(or some other interpolation rule, like linear interpolation). Then the main question 
consists in verifying if such function Uh converges somehow as /i — )• to a solution u 



(which we hope to be a classical, strong, weak, or other type of solution) to problem (5.1 ). 
To this regard, we first choose a suitable discretization in time of the source term 

/, such that the piecewise constant interpolation of this sequence produces a function 
/C^^i) (defined by means of (|5.3[)) verifies the property 



By means of these discrete approximate solutions we introduce the following notion of 
mild solution: 



Definition 5.2 We say that u G C((0,T);X) is a mild solution to (5.1) if it is obtained 
as uniform limit of the approximate solutions Uh, as h ^ 0. The initial data are taken in 
the sense that u{t) is continuous in t = and u{t) — )• uq os t — )■ 0. Besides, we say that 
u G C((0, oo);X) is a mild solution to ( |5.1[ ) in [0,oo) if u is a mild solution to the same 
problem in any compact subinterval I C [0, oo). 

In order to state a positive existence result, we need to restrict the class of operators 
according to the following definitions. 

Definition 5.3 Let A : D{A) C X ^ X be a nonlinear, possibly unbounded operator. 
Let Rx{A) be the range of I + XA, a subset of X. 

(i) The operator A is said accretive if for all X > the map I + is one-to-one onto 
R\{A) C X , and the resolvent operator Jx '■ Rx{A) X is a (non-strict) contraction in 
the X-norm (i. e., a Lipschitz map with Lipschitz norm 1). 



(ii) We say that A satisfies the rank condition if Rx{A) D D{A) for all X > 0. In 
particular, the rank condition is satisfied if Rx{A) = X for all X > 0; in this case, if A 
is accretive, we say that A is m-accretive. 
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We are now ready to state the desired semigroup generation result, that generahzes the 
classical result of Hille-Yosida (valid in Hilbert spaces and for linear A) and the variant by 
Lumer and Phillips (valid in Banach spaces, still for linear A), and provides the existence 



and uniqueness of mild solutions for problems of the type (5.1) in the case / = 0: 



Theorem 5.1 (Crandall-Liggett) Suppose that A is an accretive operator satisfying 
the rank condition. Then for all data uq E D[A) the limit 

(5.4) St{A)u^ = lim (Ji/„(^))%o. 

exists uniformly with respect to t, on compact subset of [0,oo), and u{t) = St{A)uQ E 
C([0, oo) : X). Moreover, the family of operators {St{A)}^yQ is a strongly continuous 
semigroup of contractions on D[A) C X . 

Using a popular notation in the linear framework, we could write St{A)uQ = e~*'^uo, and 



because of this analogy formula (5.4) is called the Crandall-Liggett exponential formula 
for the nonlinear semigroup generated by —A. The problem with this very general and 
useful result is that the X- valued function u{t) = St{A)uo solves the equation only in 
a mild sense, that is not necessarily a strong solution or a weak solution. Though it is 
known that strong solutions are automatically mild, the correspondence between mild and 
weak solutions is not always clear. For the FPME this issue has been discussed in detail 
in [221 [23]. 

In addition, the Crandall-Liggett Theorem result can be extended when we consider 
nontrivial source term /, according to the following result 



Theorem 5.2 Suppose that A is m-accretive. If f ^ L^{0,oo; X) and uq E D{A). Then 



the abstract problem (5.1) has a unique mild solution u, obtained as limit of the discrete 



approximate solution Uh by ITD scheme described above, as h ^ 0: 

u{t) := lim Uh{t) , 

and the limit is uniform in compact subsets of [0,oo). Moreover, u E C([0, oo);X) and 
for any couple of solutions ui, U2 corresponding to source terms f\, f^ we have 

\\ux{t) - u^mx < \\ui{s) - U2{s)\\x + r ||/i(r) - /2(r)||xdr 



for all < s < t. 

There is a wide literature on these topics, starting with the seminal paper by Crandall and 
Liggett [19], see also [18j and the general reference These notes are based on Chapter 10 
of the book [IH], cf. the references therein. The last formula we have mentioned introduces 
the correct concept of uniqueness for the constructed class of solutions. Characterizing 
the uniqueness of different concepts of solution is a difficult topic already discussed by 
Benilan in his thesis jlOj . 
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5.2 Parabolic Symmetrization 



In order to apply this theory we have to check that the operator associated to our evolution 
problem A is m-accretive or that it is accretive and the rank condition holds, in the sense 



of definition 5.3 A main question in this approach to nonlinear evolution is the corrected 
identification of the operator. This has been done in [22] as follows. 

If Mo G Li(M^)nL°°(M^), we introduce the nonlinear operator Aq : D{Ao) C L^{R'^) 
Li(M^), defined by 

Ao{u) := {-Ar/^A{u) , 

with domain 

D{Ao) := {v G L\R^) n L°^(M^) : Ao{v) G L^R^) n L^{R^)} . 



Returning to the results of Subsection |3j 1 , we see that the contractive property (3.5) 



implies that this operator is accretive in the space X = L (R ). On the other hand, 



Theorem 3.1 and its extension on R^ gives the rank condition in n L°°. By closing 
this operator with respect to the norm ||.||i we find an operator A that is m- accretive in 



Therefore, we can use Theorem 5.2 that implies that there is a unique mild solution to 
(5.1), obtained as a limit of discrete approximate solutions by the ITD scheme. In the 
case A{u) = vJ^ the extra regularity of these solutions is discussed in detail in the papers 
[221 [23]. For general ^ see [25j. 

We can now use this method to prove a symmetrization result for Fractional Fast Diffusion 
Equations, including in particular the well-known linear fractional heat equation, 

ut + {-^y/^u = 0. 



Theorem 5.3 Let u be the mild nonnegative solution of the FPME (1.5) with < a < 2, 
posed in Q = R^ , with initial data uq G L^{R^) > and nonlinearity A{u) given by a 
concave function with A{Q) = and A' [u) > for all u > 0. Let v be the solution of the 
corresponding symmetrized problem 



(5.5) 



( vt + {-Ay/^A{v) = xGR^ , t>0, 
\ v{x,0) = uf{x) xGR^ . 

Then we have for all t > 

(5.6) u*{\x\,t) '<v{\x\,t). 

In particular, we have \\u{-,t)\\p < \\v{-,t)\\p for every t > and every p G [l,oo]. 
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Proof. According to what explained before, we use the imphcit time discretization scheme. 
For each time T > 0, we divide the time interval [0,T] in n subintervals where 
tk = kh and h = T/n. We construct then the function Uh which is piecewise constant in 
each interval {tk-i,tk], by 



Uh,i{x) if t e [0,ti] 
Uh,2{x) if t G {ti,t2] 



(x) if t £ {tn-1 

where Uh,k solves the equation 

(5.7) h{-Ay/^A{uh,k) + uh,k = uh,k~i 



with the initial value Uhfi = ^o- Similarly, concerning the symmetrized problem (5.5), we 
define the piecewise constant function by 



' Vh,i{x) ift G [0,ti] 
Vh,2{x) ift G itl,t2] 



Vh{x,t) = < 
where Vh,k{x) solves the equation 

(5.8) h{-/\Y/^A{vh,k) + vh,k = vh,k-i 

with the initial value Vhfl = Uq . Our aim is now to compare the solution Uh^k to (5.7) 



with the solution (5.8). We proceed by induction. Using Theorem 3.7, we get 

If we suppose by induction that Uj^k-i ^ Vh,k-i and call Uh,k the (radially decreasing) 
solution to the equation 

h{-AYl'^A{uh,k) + Uh,k = uf^k-v 



Theorem 3.7 and Theorem |3.8| imply 

(5.9) Mufk) ^ Muh,k) -< A{vh, 



hence (5.9) holds for all /c = 1, . . . , n. Therefore, by the definition of Uh and Vh, we find 
(5.10) Aiuhi-,t)*) <Aivhi-,t)) 
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for all times t. Using Lemma 2.1 with the choice $ = F o B, where F > is convex and 
F(0) = 0, we obtain 

/ F{u*{x,t))dx< [ F{vh{x,t))dx, 

which in turn yields 

(5.11) u*i;t)^Vh{;t). 



Now Crandall-Liggett Theorem (5.2) implies 



Uh ^ u, Vh ^ V uniformly. 



Then passing to the limit in (5.11) we get the result. □ 



5.3 Symmetrization for the equation with a left-hand side 

We now consider the case / G L^{Q), Q = x (0, oo) and / ^ 0. In that case the 
semigroup generation Theorem |5.2| can still be applied to obtain the so-called unique mild 
solution of the evolution problem 



ut + {-Ay/^A{u) = f x£R^,t>0., 
u{x,0) = uo{x) xGM^. 



(5.12) 

As explained in Subsection |5.1[ we need to perform a discretization of / adapted to the 
time mesh = kh that we have used above, let us call it {fl^^}, so that the piecewise 
constant (or linear in time) interpolation of this sequence produces a function f^^\x,t) 
such that 11/ - 1 — )■ as /i — 0. Then we use the previous implicit discretization 
scheme, now in the form 

(5.13) ^{uh,k - Uh,k-i) + i-Ar/^Aiuh,k) = , 

to produce the semi-discrete function {uh{x,tk) = Uh.k{x) ■ k = 0,1,---}, which after 
interpolation in time serves as /i-approximation to the mild solution u{x, t). According to 



(5.13) we have to solve the elliptic problems 



(5.14) h{-Ar/^A{uh,k) + Uh,k = Uh,k-i + h fi'''> , 

and we can use the theory developed in Section [3j Then we have the following result 



Theorem 5.4 Let u be the nonnegative mild solution of the FFME (1.5) with < a < 2, 



posed in Q = M^, with initial data uq G L^(M^), uq > 0, right-hand side f G L^{Q), 
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/ > 0, and nonlinearity A(u) given by a concave function with A{0) = and A'(u) > 
for all u > 0. Let v be the solution of the symmetrized problem 



(5.15) 



vt + {-Ay/^A{v) = f*{\x\,t) xeR^, t>0, 
v{x,0) = uf{x) xGM^, 



where /*(|x|,t) means symmetrization of f{x,t) w.r. to x for a.e. time t > 0. Then, for 
all t > we have 

(5.16) u*{\x\,t) '<v{\x\,t). 

In particular, we have \\u{-,t)\\p < \\v{-,t)\\p for every t > and every p £ [l,oo]. 



The proof follows the lines of Theorem 5.3, so we leave the details to the reader. 

Remark. As an easy extension, we can have also a result about comparison of concen- 
trations for the radial solutions of two evolution problems, if we assume that the initial 
data satisfy the condition uq^i -< mo,2 and the right-hand sides satisfy fi{-,t) -< /2(-,i) for 
almost all t > 0. The conclusion is that ui{-,t) -< U2{-,t) for all t > 0. Let us remind the 
reader that the result holds only if A is linear or concave, as assumed above. 



6 Negative result about concentration comparison for the 
Fractional PME 

As in the elliptic case, it came to us as a surprise that the comparison result could not 
be proved for general nonlinearities A without the assumption of concavity. It turns out 
that for convex powers it does not hold. Here we will state and prove the negative result 
about concentration comparison for solutions of the Fractional Diffusion Equation in the 
range of exponents m > 1, usually known as Slow Diffusion. We first argue in a formal 
way, since we give later the justification of some details. 

• Let us consider the Fractional PME: ut + (— A)°'/^u'" = in and consider first the 
Barenblatt solution that was studied in [ST]. Here we suppose that 

m > {N - a)/N =: rrtc. 

Let us fix the mass to 1 for simplicity. The Barenblatt solution has the form 

c/i(x,t) = r"F„,i(e), i = \x\t-P 

with a = N/3 and /3 = l/(A^(m — 1) -|- o") and -Fm,i is the Barenblatt profile of mass 1. It 
is also known that as — )• oo we have 

(6.1) Fm,i{0-cr^''+'^^- 
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This means that for large x and t ~ (so that ^ ~ oo) we get 

Ui{x,t) r^Ct^\x\-^^+'^\ A = /3a. 

This approximation holds uniformly for all \x\ > C large and all < t < r if r is small 
enough, that is the error is higher order small 

Ui{x,t) = Ct^\x\-^^+''\l + e). 

Now, let us choose the initial data to associate to the equation. Let 

«o(x) = 0(|x|)>O, 

where (/> is a smooth and compactly supported function in the ball of radius one, having 
mass 1. Suppose also that is rearranged. Let tis call u the solution to the equation with 
such choice of the data. We have that (— A)'^/^Uo' is a bounded function that decreases 
at infinity like Ci By virtue of the equation we have 

ut(x,0) = -(-Ar/2^,- 
so that for small t we have approximately 

u{x,t) ~ -ti-^f/^u"^ 

and this behaves as t — )• oo like Ct\x\~^^~^'^\ A rigorous proof of this behaviour will be 
given in the next two subsections. 

The conclusion is that the concentration comparison result is not true. The reason is 
that the exponent 

^ = ^^= a + Nlm-1) 

is less than 1 (precisely for m > 1). Clearly, Uq = S{x) is more concentrated than uq, but 

u{x,t) is larger than U{x,t) for very large a: if t is quite small. This is incompatible with 
satisfying the concentration comparison and having the same mass. 

In the graphics of Figure 1 we show the relative evolution of two solutions in time. 
Initially the concentrations are ordered, later they are not. The parameter m is 2. 

6.1 Supersolution. First tail estimate 

Let us now give a rigorous derivation of the tail behaviour. Fist, we have a preparatory 
step. 

Theorem 6.1 Let u{x,t) he a classical solution of the FPME with initial data, uq{x) > 
such that uo{x) < 1 in the ball -Bi(O) and uo{x) < |a;|~^^+'^^ for \x\ > 1. Then there is a 
time ti > such that 

(6.2) u{x,t) < 21x1-^^+"^ 

if \x\ > R and R is large enough and <t <ti. 
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Figure 1: Comparison of FPME evolution at four consecutive times. 



Proof. We consider the FPME for m > 1 and initial data uq is 1 in the ball of radius 2. 
Then for all times the solution will be bounded by 1. 

We want to construct a super-solution of the form 
(6.3) Uix,t) = il + bt)F{\x\), 

where F > has to chosen. I will need F{r) ~ Cir~^^~^'^^ as r = |x| — )• oo to get the 
desired conclusion after the comparison argument in the following way: 

u{x,t) < U{x,t) < 2F(r) < 2Cir~(^+'") 

if r is large and t < ti/b 0. 

To establish such comparison we first note that Ut = bF{x). Using the notation 
, we also have 

L,[/™ = (1 + bt)L,F'^{x) 

As we have pointed out, F{r) ~ Cir'(^+'^) as r oo so that L^F"' = 0(r-(^+'^)) for 
r > 1 (cf. Lemma 2.1 of [E])- It follows that there is a constant /c > such that 

kF + L^F"^ > 0, everywhere in M^. 

Therefore, we will have 

Ut + L^U"" = 6F + (1 + btri^F"" > 
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if 6 > k{l + bt), i.e. if 6 > and t < (b — k)/kb, for instance for b = 2k and t < 1/6 = l/2k. 

Under such assumptions, the viscosity method will work in the exterior region Q = 
{{x,t) : > 1,0 < t < ti}, and this will prove that U{x,t) < u{x,t) in Q as desired. 

We finally check the application of the viscosity method. Indeed, the boundary condition 
at r = 1 is 

U{l,t) > F{1) > 1. 

so U{x,t) > u{x,t) on the lateral boundary of Q located at r = 1. Same comparison 
is trivial for t = 0. We only need to argue by contradiction at the first point where 
the classical solution u touches U from below to conclude that u{x, t) is strictly less 
than U{x,t) in Q. The contradiction at the point of contact is explained in [13]. The 
construction of classical solutions is done in 



Remarks, (i) We have done the argument for R = 1. If we want to change the radius to 
> 1 we may use the scaling of the equation. 



(ii) Lower estimates that match the tail behaviour (6.2) are derived and used in [38j 



6.2 Supersolution. Sharp tail estimate 

Theorem 6.2 Let u{x,t) be a classical solution of the FPME with initial data uq{x) > 
such that uo{x) < 1 in the ball -Bi(O) and uo{x) = for \x\ > 1. Then there is a time 
ti > and constants C* and R such that 

(6.4) u{x,t)<C*t\x\-^^+'''^ 

if \x\ > R and R is large enough and < t < ti. 



Proof. We proof is a delicate variation of the preceding one. We still consider the FPME 
for m > 1 and initial data uq is 1 in the ball of radius 2. For all times the solution will be 
bounded by 1. We consider a supersolution of the form 

C/"(x, t) = G{x) + F{xr, 

where i*" > is chosen as before. I take F{r) ~ Cir^^^"*""^^ as r — > oo. Again, it follows 
that there is a constant k > such that 

kF + L^F'^> 0, everywhere in M^. 

Next we choose G > compactly supported and such that Lq-G = cq > on the support. 

As r — )• oo, we get the usual L^jG ~ _ cr^('^+^^). 

We get the formula 
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which reduces to Ut = bF when G = 0. In any case it is nonnegative, Ut > 0. We also 
have 

Then when G = we will have 

Ut + L<,C/™ = bF + (6t)'"L<,F™ + L^G > 

if b > Bq, b > 2k{bt)"^, i.e. if b"^~^t'^k < 1/2, which imposes a condition above on t. 
On the other hand for G > we have 

Ut + L^U"" > CO + (UrL^F^ > 

if G2{bty"- < cq. Both conditions are fulfilled if < t < ti. 

Is this is the case the viscosity method will work in the region Q = > 1,0 < t < ti}. 
and this will prove that U{x, t) < u{x, t) in Q. 

Indeed, the boundary condition at r = 1 is 

U{l,t) > G{1) > 1. 



Remarks. The rate of decay (6.4) of the tail of such solutions at infinity is optimal as a 
consequence of the construction of suitable sub-solutions with the same exponents in the 
X and t dependence, which is done in [38] . 

The counterexample is heavily technical. Surprisingly, the situation becomes much clearer 
when we let m ^ oo. This is studied in 1521. 



7 Comments, extensions and open problems 

-In a companion paper [53j we will use symmetrization results of this paper to obtain 
sharp a priori estimates with best constants for some functional embeddings involving the 
solutions of the linear fractional heat equation or its fast diffusion relative. 

-As an extension of the above results, we could consider equations that involve a more 
general version of the fractional Laplacian operator, in the same way that the standard 
symmetrization applies to elliptic equations with coefficients. 

-The elliptic and parabolic counterexamples have been constructed for the problems posed 
on the whole space. They could also be constructed for solutions defined on a bounded 
domain, say a ball, with zero Dirichlet boundary conditions. The argument is as follows: 
we consider the problems posed in a sequence of balls Bji expanding so that i? — t- oo with 
same data of compact support. According to [22l [23] the solutions ur converge to the 
solutions of the limit problem in the whole space. Now, for the limit equation there is a 
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counterexample. We deduce that there is a counterexample before the limit. We leave 
the details to the reader. 

-Another interesting problem would be obtaining a priori estimates for solutions of elliptic 

and parabolic problems of this type with Neumann boundary conditions using symmetriza- 
tion techniques. A good indication is that conservation of mass is true for both elliptic 
and parabolic problems. 

Let us now list some open problems that have arisen in the course of the work: 

-We do not know how to deal with concave nonlinearities A in bounded domains. 

-We do not know how to do the elliptic or parabolic comparison in the case of more 
general function A, if it is neither concave or convex. 

-Finally, we wonder if there is a partial or alternative theory that replaces the failure of 
the concentration comparison result for the fractional porous medium equation, i.e. the 
equation 

dtu + (-A)'^/^"^ = 

with m > 1. 
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